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Abstract 

We study the dynamics of an electron confined in a one-dimensional double quantum dot in the 
presence of driving external magnetic fields. The orbital motion of the electron is coupled to the spin 
dynamics by spin orbit interaction of the Dresselhaus type. We derive an effective time-dependent 
Hamiltonian model for the orbital motion of the electron and obtain a synchronization condition 
between the orbital and the spin dynamics. From this model we deduce an analytical expression 
for the Arnold tongue and propose an experimental scheme for realizing the synchronization of the 
orbital and spin dynamics. 
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I. INTRODUCTION 



Phase synchronization and related phenomena are among the most fascinating effects of 
nonhnear dynamics. Besides the deep fundamental interest [l|-l7| , phase synchronization has 
a broad range of applications in chemistry js | , ecology {q], astronomy 10 1, in the field of 



information transfer using chaotic signals ll|, and for the control of high frequency electronic 
devices 12| . In nonlinear dissipative systems, phase synchronization occurs if the frequency 



of the external driving field is close to the eigenfrequency of the system. In this case, for 
a certain frequency interval of the driving field the oscillations of the nonlinear dissipative 
system can be synchronized with the perturbing force. Usually, for stronger driving fields 
the frequency interval for the synchronization becomes broader and the synchronization 
protocol is more efficient. The broad and growing interest in the phase synchronization calls 
for the analysis of new possible realizations of this phenomenon. A particularly interesting 
issue is the application of the synchronization protocols for magnetic nanostructures, which 



have rich applications ll-[l6| and exhibit interesting nonlinear dynamical properties that 
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18| 



can be exploited as a testing ground for dynamical systems 

A principal challenge in nanoscience is to find an efficient procedure for the manipulation 
of the systems states. A high level of accuracy on the state control is required especially in 
such applications as quantum computing, where a precise tailoring of the entangled states 
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20j. With this in mind several physical systems were considere d up 



is highly desirable 

to now, e. g. Josephson junction qubits and Rydberg atoms in a quantum cavities 



ion traps 
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2l|, single molecular nanomagnets 22|, and nanoelectromechanical resonators 



24j . Among others, one of the most promising systems are electron spins confined 



in two-dimensional quantum dots [25 



based quantum dots 



27| and in one-dimensional nanowires and nanowire- 



28l-l3l|. The key element of the corresponding models is the spin orbit 



(SO) coupling term, which is linear in the electron momentum. Such momentum-dependent 
coupling offers a new way of manipulating the spin by changing the electron momentum via 
a periodic electric field. This is the idea of the electric-dipole spin resonance propose d by 
Rashba and Efros for the electrons confined in nanostructures on the scale of 10 nm 27 |. 
However, the external electric field can strongly affect the orbital dynamics and thus the 
system is driven out of the linear regime 32| . The nonlinearities usually result in a complex 



behavior of the affected systems and their dynamics might become complicated and even 



unpredictable. On the other hand, the nonhnearity may also lead to a number of interesting 
phenomena. In spite of the huge interest in the systems with SO coupling, the influence 
of the spin dynamics on the orbital motion was not addressed yet in full detail. With the 
present work we would like to bridge this gap. Our goal is to investigate the possibilities of 
controlling the orbital motion of the electron via external magnetic fields acting on its spin 
and via the spin-orbit coupling. That can be considered as opposite to the electric-dipole 
spin resonance protocol proposed by Rashba and Efros We will demonstrate that: (i) 
by using an external driving field one can achieve a sufficient degree of control over the 
orbital motion, and (ii) as a result, one can devise a very efficient synchronization protocol 
of the orbital motion and the spin dynamics based on the application of a pulsed external 
magnetic field. 



II. THEORETICAL MODEL 

We consider a model system of a single electron confined in a double quantum dot de- 
scribed by a potential of the form U{x) = Uq [—2 {x/d)^ + (a;/(i)^]. Here Uq is the energy 
barrier separating two minima with 2d being the distance between them. We assume the 
system is dissipative, and the dissipation is a thermal effect appearing due to a coupling 
to environment. The dissipation, which impacts mainly on the orbital motion, is essen- 
tial for the synchronization processes we are going to discuss later in the text. For strong 
driving magnetic fields, the influence of the dissipation on the spin dynamics is negligibly 
small and can be ignored. In addition, we assume that the temperature is low enough to 
prevent the activated over-the-barrier motion. For the particular value of Uq ~ 20 meV the 
low-temperature regime means T < 100 K. For the GaAs-based structure with the electron 
effective mass m being 0.067 of the free electron mass and d ~ 100 nm, the tunneling prob- 
ability is small and a classical consideration is justified. To quantify the SO interaction we 
use a coupling term of the Dresselhaus type Hso = otPxCy^-, where Px is the momentum of 
the electron and is the Pauli matrix. Therefore, the Hamiltonian of the one dimensional 
system reads: 

H = ^ + U{x) + aP^a" + fiBgB.it)^ + I^BgBxit)^, (1) 
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FIG. 1. Schematic illustration of the infinite series of external magnetic field pulses applied to the 

oo 

system. A series of the short pulses Bz{t) = Bq 6r{t — nT) with the pulse width r and time 

i=0 

interval between pulses T, is applied along the z axis. Series of the pulses with the larger width T 

oo 

and a shorter interval between the pulses r, Bx{t) = BqYI ^r(^ ~ tiT') is applied along the x axis. 

t=o 

The amplitude of the pulses Bq is the same in both cases. 



where /is is the Bohr magneton and g is the electron Lande factor. Here Bz{t) = Bq"^ 6r{t— 

t=o 

nT) is an infinite series of external magnetic field pulses with the pulse strength Bq, which 
is applied along the z axis. The temporal width of the pulses applied along the 2;-axis is 
smaller than the interval between the pulses r <^ T (in what follows we set T = 1). On the 

oo 

other hand, for the pulses along the x-axis Bx{t) = Bq^ Ar(t — rn) the pulse duration is 

t=o 

larger than interval between pulses T <^ t (cf. Fig. 1). A different route to the control of 



the spin-djTiamics in double quantum dots via electric field pulses is outlined in 
Introducing the characteristic maximum momentum of the electron p™'*^ 



33| 



we can estimate the maximal precession rate of the spin due to the SO coupling Ql 



)max 
^so 



{2a /h) y/2mUo, while the magnetic field pulse of the amplitude B = Bq induces a spin 
precession with the rate Qb = f^BldlBo/h. Therefore, if Qb > ^^'^^ we can during the pulse 
neglect the spin rotation produced by the SO coupling and the spin is completely controlled 
by the external driving fields. We need a protocol with two driving fields in order to fulfill 
the synchronization requirements as discussed later in the text. Namely, for the control of 
the spin dynamics via the external driving fields, the amplitudes of the fields should be large. 



Bq > (2a//iB|fi'|)A/2mf/o. On the other hand, a strong constant magnetic field produces a 
high frequency precession of the spin Qb = f^B\g\Bo/h, while for the synchronization we 
need to tune the precession frequency up or down keeping fixed the strong driving field 



amplitude. Below we will show that the optimal conditions for the synchronization are 
realized using two types of the driving pulses. Applying short pulses along the z-axis r ^ T, 

oo oo 

B;^(t) = Bq"^ 6r{t — riT) and long pulses along the x-axis B^it) = BqYI ^rii — tu), we 

t=Q t=o 
can realize a spin precession with a frequency, which is inversely proportional to the time 

interval between the short pulses f2 ~ 1/T independently from the driving field strength 

Bq. In what follows, for convenience we use dimensionless units via the transformations 

E E/AUo, X x/d, t t^AUo/m, P^/ ^fhnU'Q, e a/AUo. 

III. DISSIPATIVE SYSTEM AND THE PROBLEM OF PHASE SYNCHRONIZA- 
TION BETWEEN ORBITAL AND SPIN MOTION 

A. Spin dynamics in pulsed magnetic fields 

As was stated above the synchronization can occur if the frequency of the driving field 
is close to the eigenfrequency of the nonlinear dissipative system. If this is the case, in the 
particular frequency interval, the oscillations of the nonlinear dissipative system and the 
field can be synchronized. With the increase in the driving field amplitude, the synchroniza- 
tion can occur in a broader frequency interval, and the synchronization protocol becomes 
more efficient. Our aim is to develop a method for the synchronization of the dynamics 
of the electron spin and the orbital motion, using an external driving magnetic field and 
SO coupling. Although the magnetic field is not coupled to the orbital motion directly, a 
sufficiently strong field infiuences the orbital motion through the spin dynamics if the SO 
coupling is present. If SO term is relatively small fl^^^ < VIb, that is 

ns = ^^^>n-='-^V^o, (2) 

the spin and, correspondingly, the orbital motion can be controlled externally. From Eq. ([T]) 
it is easy to see, that in between the short pulses the electron spin rotates around the x-axis 
and the equations of motion for the electron spin in this case read 

= 0, &y = -VLBa\ a' = fi^a^. (3) 

On the other hand, during the short pulses we have 

= -nea^, &y = VtBa\ = 0. (4) 
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Considering the dynamics due to the pulse acting on the spin at the moment of time t = to, 
we can spht the evolution operator Tcv defined as 



a(t'-^+T)=n^cT[tl^) (5) 
into two parts: T^v = x T^, where 

a{t^^^+T)=fnrT{t^:^). (6) 

Here we introduced the notations t\^^ = to + and = to — The operator Tr describes 
the rotation of the electron spin around the x-axis produced by the long pulse of the external 

oo 

driving magnetic field B,j.{t) = ^rit — tu), which is applied along the x -axis and T^. 

t=o 

oo 

corresponds to the evolution produced by the short pulses Bz{t) = Bq^ 6r{t — nT) applied 

t=o 

along the ^-axis. Integrating Eq. (|4]) for a short time interval t e (^^cT^^^o^'j we obtain 

fsia^ = a%t[-^) cosiQBT) - ^^(t!,-') sin(fiBr), (7) 
fsia^) = a^'(4"^) sin(fiBr) + ^'^(t'^) cosiQBT). (8) 

Integrating Eq. (|3]) during the time interval t G ( ^o^^, ' + T ) of long applied pulse we find 



T«K) = Wl - a- ti,+M cos{QbT), (9) 



Tn{a') = a/1 - [a- C sm{nBT). (10) 



Combining Eq. ([7]) with Eq. (|9]) we finally can reconstruct the complete picture of the full 
time evolution of the electron spin: 



^n+l = v/l-K+l)'cOs((n + l)nBT), 

<+i = v/l-«+iPsin((n + l)fi^T), (11) 



'n+l 



(T^cos(fiBi") — a/I — {a'nY sm.{VtBT) cos{nVtBT). 



The recurrent relations Eq. fITT]) describe the spin dynamics. Accuracy of the employed 
approximations may be checked by the validity of the normalization condition cr^ = 1. In 
order to identify, whether the nonlinear map ( ITT]) is chaotic or regular, we evaluate the 



Lyapunov exponent for the spin system 



34| . Taking into account the peculiarity of the 




FIG. 2. Lyapunov exponent as a function of the initial spin component ctq- After A'^ = 1000 
iterations the Lyapunov exponent is negative A(o"o ) < 0, meaning that the spin dynamics is regular. 
T = 1, Qbt = I, 0,bT = 20. At these conditions, the dependence is weak. 

system (fTT]) . which is the fact, that the equation for the x component dn+i is self-consistent 
o'n+i = fi'^n)y deduce for the Lyapunov exponent 



hm — In 



Sa 



hm — In 



(12) 



Here, cTq is the initial value of the spin projection and the small increment of the initial 
values (Tq + Sa quantifies the sensitivity of the recurrence relations ffTT]) with respect to the 
slight change in the initial conditions. After some algebra from Eqs. ffTTj) and f|T2|) we finally 
obtain: 



N-l 



lim — In 

N^oo N ^ 



n=0 



da't 



N-l 



N ^ 



n=0 



cos(r2 



B 



: sin(f2B) cos(nr2By 



(13) 

The results of the numerical calculations are presented on Figs. [2] and [31 From Fig. |2] we 
see that the Lyapunov exponent is negative, A((Tq) < and therefore the spin dynamics 
is regular, since the initial distance between two neighboring trajectories starting from the 
initial points and ctq + ba is not increasing asymptotically after an infinite number of 
iterations ^cre^'^'-'^''^ Therefore, from Figs. [2] and [3] we conclude, that the dynamics of 
the electron spin is controlled by the magnetic field pulses, thus following equation o'-^if) = 
(Tq cos (fit) we achieve the spin manipulation by magnetic fields. The spin rotation frequency 
is determined by the time interval between the short pulses Q ^ 1/T. 
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FIG. 3. Time dependence of the spin projection cx^ = (T^(n). T = 1, 0,bt = 1, ^IbT = 20. From 
Fig. Owe see that = cr^(t) is periodic in time and the frequency of oscihations il. is inversely 
proportional to the time lapse between pulses Q ~ 1/T. 

B. Synchronization of the spin and the orbital motion. 

With the spin dynamics discussed in the previous subsection, for the orbital motion of 
electron we can write the following effective Hamiltonian assuming that cr^'(t) ~ cos{Qt): 

p2 

H=^ + U(x) + aP^ cos(fit). (14) 
2m 



The equation of motion for the system f[Hj) has the form: 

X + ■yx — X + x"^ = —(3 sin(f2i:) , (15) 

where two new dimensionless quantities are introduced: (3 = aQm/AUo and 7 — )■ 'y/AUom. 
We seek a solution of Eq. (ITSil using the following ansatz 

x{t) = lA{t)e'^' + lA*{t)e-'^' . (16) 

Assuming that the amplitude A{t) in Eq. ( fT6l) is a slow variable the following condition 
applies 

i(t)e^™ + i*(t)e-^™ = 0. (17) 



1 2 3 4 5 6 

Y 

FIG. 4. Arnold tongue diagram in terms of the decay constant 7 and field frequency fi, plotted 
using the conditions ()25p . Shadowed regions define parameter values for which the synchronization 
is possible. 

Taking into account Eqs. ( |T6l) and ( |T7j) from Eq. (fT5l) we deduce: 
zfii(t)e'™ - A(t)e*™ + c.c.^ + 

+7 (^y A(t)e*^* + c.c.^ - QA(t)e*^* + +c.c.^ + (18) 

+- (A3(t)e3'™ + 3U(t)PAme*™ + c.c.) = ^ (e^™ - c.c.) . 
8 ^ 2i ^ ^ 

Multiplying Eq. f ll8p by the exponent e"**^* and averaging it over the fast phases we find: 

Introducing the notations 

m = 2^z(t), r=^. A = — . . = ^. (20) 

we can rewrite Eq. flTIJl) in a more compact form 

i(r) + zA^(r) = -z{t) + — k(r)pz(r) + e. (21) 
Inserting 2;(r) = i?(r)e*'^*''^-* into Eq. fl2T]) . for the real and imaginary parts we obtain 

R{t) = —R + e cos(p{t), 

0(r)+A = ^i?2(r)--^sinv,(r). ^^^^ 



^ 3- 
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FIG. 5. Arnold tongue in terms of the parameters A = (^2^ + l) /'^VL, (A, e) plotted using Eq. ()26p . 
Point h belongs to the domain where a synchronization is not possible, while a belongs to the 
synchronization domain. 

Using Eq. (!22|) and setting i? = 0, = for the stationary solutions we obtain: 

3 



/(0 = e + e(^^e-Aj =e\ (23) 

6,2 = f (2A ± v/A^) , (24) 

where = ^ and are roots of the equation df{^)/d^ = 0. In order to identify the Arnold 
tongue 351], which shows the regions in which a synchronization is possible, we utilize the 
standard condition df{^)/d^ = 0. From Eq. ( 12^ it is not difficult to see that the roots of 
the equation df{^)/d^ = are real if A > v^. Taking into account that A = (il"^ + 1) /'yO, 
we can rewrite the inequality in the form fi^ + 1 > 7f2-\/3- Consequently we obtain the 
following criteria for the synchronization 

0<fi<i(7V^- V372-4), ^(7^3+ a/372 -4), 7 > (25) 

Graphical representation of the conditions ( l25l) is shown in Fig. 4. 

Eq. (25) defines the synchronization area in terms of the external field frequency fl 
and the dissipation constant 7. The minima points of the function /(O, (23) does not 
depend on the SO coupling constant a. Therefore criteria (25) is independent of the values 
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FIG. 6. Orbital and spin dynamics plotted using the exact numerical integration of the Eq. (jlSp . 
Panel a) corresponds to the values of parameters for which a synchronization is possible (see Fig. [5l 
point a). In particular ^e/\/0 = 4.41, A = 1.88^. The right panel b) corresponds to the point b 
on FiglU ^e/v^ = 1.05, A = 3.03^, i.e. the dynamics occurs outside of the synchronization area. 
We see that in the former case the orbital and the spin dynamics are in phase and in the latter 
case the phase difference is about 7r/2 . 

of the SO coupling strength as welL Nevertheless, inserting the roots .^1,2 of the equation 
dfiO/d^ = into Eq. ( l23l) one can derive more illustrative and precise criteria in the form 
of the parametrical curve: 

2 



-(±3VA2 - 3 + 9A + A^T AVA2 - 3) 



0. 



(26) 



The parametrical curve defined by Eq. (!26l) represents the border of the synchronization 
domain, see Fig. [51 Taking into account that (3 = a^lm/ AUq, e = /3/2f27'^/2 obtain 

e am 



^ v^8f/o7=^/2- ^^^^ 
From Eq. fl271) we see that the parameters of Fig. [5] depend on the oscillation frequency 
that can be easily controlled by tuning the time interval between pulses Q ~ 1/T. All 
other parameters in Eq. flTTI) . such as the SO coupling constant a, barrier height Uq, and 
the electron effective mass m are internal characteristics of the system whereas the decay 
constant 7 is related to the thermal effects. Using Eq. ( 126|) and Fig. |5] one can synchronize 
the electron orbital motion with its spin dynamics. 
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IV. CONCLUSIONS. 



We have investigated the classical electron dynamics in a double dot potential, with the 
spin of electron being controlled by external magnetic fields. We have shown that the orbital 
electron dynamics can be controlled very effectively by the field in the presence of a spin-orbit 
coupling. Using the proposed protocol of magnetic field pulses of different duration we have 
shown that it is possible to synchronize the spin and the orbital motion of the electron. In 
particular, if the driving field amplitude is large enough, Bq > (2a;//iB|5'|)v^2mf7o, the spin 
dynamics is periodical in time. Then cr^{t) = ctq cos(r2t), where the oscillation frequency 
is inversely proportional to the time interval between pulses Q ^ 1/T and can be tuned 
independently from the amplitude of the pulses Bq. As a consequence the orbital dynamics 
can be studied with reduced effective, time-dependent, one-dimensional model. By using this 
model we found the synchronization condition between the orbital and the spin dynamics. 
Furthermore, we derived an analytical expression for the Arnold tongue that defines the 
values of the parameters for which a synchronization is possible. Since in the designed 
protocol the spin precession rate is determined by the interval between the applied pulses 
we believe that it can be realized in future experiments on semiconductor quantum dot 
devices. 
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